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■ Abstract. Each spacecraft in the Laser Interferometer Space Antenna (LISA) houses 

a proof mass which follows a geodesic through spacetime. Disturbances which change 
the proof mass position, momentum, and/or acceleration will appear in the LISA 
data stream as additive quadratic functions. These data disturbances inhibit signal 
extraction and must be removed. In this paper we discuss the identification and fitting 
of monochromatic signals in the data set in the presence of data disturbances. We also 
, present a preliminary analysis of the extent of science result limitations with respect 

' to the frequency of data disturbances. 
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H ■ 1- Introduction 

The Laser Interferometer Space Antenna (LISA) will detect and study in detail 
gravitational waves from cosmological sources [I . . Probably the most interesting sources 
that LISA will detect are binaries containing massive black holes. Merging galactic 
structures may contain intermediate mass or more massive black hole binary systems. 
As the binary inspirals it will radiate gravitational waves at low frequencies, typically 
below 10 niHz for much of the inspiral phase. Systems involving supermassive black holes 
will radiate at still lower frequencies, down into the (J.Hz frequency region. Gravitational 
wave signals such as these will help to answer questions about galaxy formation and 
massive black hole production. 

The LISA constellation will orbit the Sun some fifty million kilometers away from 
the Earth. Although far from disturbances due to the Earth-Moon system, the LISA 
environment is not completely benign. There are two environmental considerations for 
LISA that will be especially important to us: cosmic rays and large solar flares. The 
probability of a micro-meteorite encounter with one of the spacecraft is very small, and 
it is unlikely that collisions will be significant compared to the other two events. 

Cosmic rays from the galaxy and Sun bombard each LISA spacecraft continuously. 
Cosmic rays with energies above about 100 MeV will be able to travel through the 
spacecraft shielding and deposit some of their energy and charge in the enclosed proof 
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mass. Charging of the proof mass is a problem due to the resulting electrical and 
magnetic forces on the proof mass. As a gravitational wave antenna, the sensitivity 
of LISA is strongly dependent on each proof mass following a geodesic. To counteract 
charging of the proof mass, it is planned to continuously or periodically discharge the 
proof mass using a discharging lamp which shines ultraviolet light onto the proof mass 
and/or its housing, ejecting photoelectrons. 

Throughout the 22 year solar cycle there are times when the Sun is particularly 
active. At these times, the Sun will occasionally produce large solar flares. These solar 
flares produce a large increase in solar cosmic rays. During every large solar flare event, 
and possibly during discharging, the proof mass will get a momentum kick. During the 
event the proof mass will be slowly pushed in some direction and viable data may not 
be obtained. At the end of the event the proof mass will have acquired a new position, 
and will have a new velocity. 

In addition to the environmental disturbances, the telecommunication antenna 
(which communicates with Earth) periodically (perhaps as often as once per week) 
will be rotated by 0.1 radians or more. This movement temporarily will disturb the 
spacecraft and the disturbance reduction system, causing the proof masses to have 
changes in their relative positions and velocities. After the antenna has been rotated 
there will be a new gravitational field around the proof mass which will cause the proof 
mass to have a new acceleration. This is because the mass distribution of the antenna 
will not be known well enough to calculate the change in the gravitational acceleration 
accurately. 

The changes in the positions, velocities, and accelerations of the proof masses will 
appear in the LISA data stream. The new positions will change the number of laser 
light wavelengths between each proof mass and therefore will be an offset in the LISA 
data stream. A new velocity of the proof mass will appear as a slow linear ramp of the 
data, and an acceleration will appear as a quadratic term. 

These complications of the LISA data stream will hinder data analysis if not 
removed. Since we will have no other way to know the changes of the positions, 
velocities, and accelerations of the proof masses to sufficient accuracy, we must fit for 
these parameters in the data. In this paper we demonstrate a fairly simple method for 
locating a gravitational wave signal in the LISA data stream in the presence of a data 
disturbance. In £|2we describe the construction of our simulated LISA data stream; in 
J3we describe our search algorithm, and the methods used for locating the signals and 
disturbances; in 2] we apply our search algorithm to our simulated LISA data stream; 
in 33 we analyze the affect of multiple disturbances on the science results. 

2. Simulated LISA Data Stream 

Our simulated LISA data stream consists of noise, a monochromatic signal, and a 
disturbance. This section describes each of the three components in detail. 
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2.1. Noise Spectrum Realization 

We use a suggestion for the LISA sensitivity curve described in [2j which extends down 
to 3 |j.Hz. We further extend this curve to 0.3 \xHz. Our noise realization also includes 
an approximation for the binary confusion noise present in our galaxy |3] . Above 3 mHz 
the sensitivity is determined mainly by white noise in measuring the distances between 
the proof masses, and is dominated by shot noise in the LISA lasers. The spectrum 
from 3 mHz down to 0.1 mHz is dominated by the binary confusion noise. This part of 
the spectrum rises steeply at first and then roughly as 1/f from 1.6 to 0.1 mHz. Below 
0.1 mHz the spectrum is dominated by acceleration noise, with a power law exponent 
of —5/2. Thermal and proof mass charge fluctuations increase the acceleration noise 
at frequencies of 10 to 3 |J.Hz with another power of yff. It is assumed here that the 
spectrum continues to rise as l// 3 down to 0.3 |xHz. 

We create a realization of the LISA sensitivity curve in frequency space by 
generating gaussian deviates to multiply each amplitude. Inverse Fourier transforming 
the spectral amplitude gives us a times series. Our time series contains 2 15 points 
sampled at 100 second intervals, making the observation time just under 38 days. 
Throughout this paper ti is the time stamp at measurement i, and yi is the simulated 
LISA data stream with disturbance and signal at £j. 

The bulk of our analysis used one realization of the LISA data stream. To gain 
confidence in our search algorithm we produced 35 different realizations of the LISA 
data stream. If not otherwise mentioned, all results refer to the first realization. 

2.2. Signal and Disturbance 

For this analysis we have chosen to look at single monochromatic sinusoids with constant 
amplitude. The actual LISA signals may be considerably more complex. For instance, 
there may be a large number of separable signals below 100 |a.Hz, and many may change 
their frequencies considerably in one month. However, information for the single signal 
case is expected to give some indication of what is likely to happen in more complicated 
cases. 

We analyze two single signal cases. The first signal has a frequency of 3 |j.Hz, which 
is near the "limit" of the LISA sensitivity. This signal contains nearly 10 complete cycles 
within our 38 day period. In view of the low frequency, this signal should give us strong 
indications of limitations on the frequency of data disturbances. 

The second case has a signal at 100 |j.Hz. As the LISA constellation orbits the Sun 
gravitational wave signals will be modulated in amplitude and frequency. However, the 
wavelength corresponding to 100 |J.Hz is 20 AU, so the signal will have minimal spread in 
frequency space. Therefore we are justified in approximating a signal at this frequency 
as monochromatic. 

Our signals are characterized by an amplitude, a frequency, and a phase. We 
represent the signal symbolically as a function h(ti\ A, f, <ft). We define a disturbance as 
an interruption in the data stream where no data is present, and where at the end of the 
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interruption the proof mass may have a new position, velocity, and acceleration. The 
interruption has a duration of length L. The disturbance is characterized by a starting 
time, an ending time, and three parameters [p, u, a] which are the additional position, 
velocity, and acceleration of the proof mass. We represent the disturbance symbolically 
as a function g(ti]p,u,a). 

Figure ^ shows one realization of the LISA data stream with a monochromatic 
signal at 3 |a.Hz and a disturbance located at 11.6 days after the beginning of the data 
set. 



2 




-1 t , , 

10 20 

sample time (days) 




-5 -4 
Log 10 Frequency (Hz) 



Figure 1. The top frame depicts one realization of the LISA noise spectrum in 
the time domain. The center frame adds to the noise a monochromatic sinusoid with 
(A, f, </>) = (2x 1(T 17 , 3uHz, 0°), and a disturbance with [p, u, a] = [lCT 17 , KT 23 , 1(T 28 ] 
located at 11.6 days with an interruption length of 17 minutes. The bottom frame 
shows the spectral amplitudes of the top two frames, the dotted line is our LISA 
sensitivity curve. 



3. Search Algorithm 

We use an iterative approach to finding the signal and disturbance parameters. Our 
algorithm utilizes both the downhill simplex method of Nelder and Mead Ej and 
the directional downhill decent method of Marquardt [6 . Given an initial collection of 
possible solutions to the search problem (generated randomly) we use the Nelder and 
Mead method to quickly find a local minimum of the problem. 
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As explained in [3] the Nelder and Mead method is very useful for producing a 
rapid initial decrease in estimator values in low parameter dimensions. But in higher 
dimensions the Nelder and Mead method fails to find minimized estimator values. 
However, the Nelder and Mead method can be forced to continue minimizing the 
estimator with prodding. Some of our prodding is randomly restarting our simplex 
centered about the previous solution. Other prodding involves modifying the estimator. 

Once the Nelder and Mead method ceases to refine the solution we turn to 
the directional downhill decent method of Marquardt. This search method uses 
derivative information to minimize the estimator. Finishing with the Marquardt method 
guarantees a minimized estimator, typically with a value much less than that provided 
by the Nelder and Mead method. 

3.1. Evaluation Criteria 

The search algorithm is dependent on an evaluation criterion which determines how 
likely a solution is to be the answer. Sine we have correlated noise, we cannot determine 
an uncertainty associated with each measurement point in the time domain. But we 
do have some knowledge of the noise: we know its frequency spectrum. Therefore 
we can construct a \ 2 in the frequency domain with a diagonal correlation matrix: 
X 2 = Eilo 1 [ r (/i)/°"(/i)] 2 > where r(/;) is the Fourier transform of yi - h(U) - g(U), and 
cr(fi) is our LISA sensitivity curve. 

There are two modifications we make to our estimator during our search. The 
first is to window the data before taking the Fourier transform. The discrete Fourier 
transform of a sinusoid is a delta function in frequency space with large wings dropping 
off as 1//. Using a modified Hanning window we suppress the wings of the function, 
but we also spread the peak of the function out to cover more frequency bins. This 
modification helps the Nelder and Mead method to locate the frequency and amplitude 
of the sinusoid in our data stream. 

The second modification we make to our estimator is to weight the residuals. After 
the Nelder and Mead has found decent parameter values we no longer window the data. 
Our search algorithm can get away with minimizing the x 2 by lowering the frequency of 
the sine wave, which has the effect of reducing the transform wings at high frequencies. 
Since there are many more points at high frequencies than there are at low frequencies 
the value of \ 2 will be reduced, but there will be extra residuals at low frequencies, since 
the sine wave frequency is not being fit appropriately. To force our algorithm to fit the 
sine wave frequency we weight the residuals at low frequencies more than the residuals 
at high frequencies. We justify this weighting by the fact that we know we are looking 
for a low frequency sine wave, so we preferentially want to reduce the \ 2 contribution 
from the low frequency residuals. 

In x 2 -tests a reduced x 2 °f 1 indicates a good fit. We have 32768-point data sets, 
with 6 parameters to be fit, which yields 32762 degrees of freedom. Any value of \ 2 
around this value is considered a good fit. But how good is good? We would like 
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some supercriterion to determine how well the search algorithm is actually performing. 
Our supercriterion compares the algorithm's fit for the signal, in the presence of both 
noise and a disturbance, to the algorithm's fit for the signal in the presence of noise 
only: A 2 = Y^[h(t; A, f, <f>) — h(t; A, f,4>) nodis } 2 . In this way we can judge how well 
different search criteria and search algorithm's remove the disturbance from different 
data streams. 

3.2. Disturbance Identification and Removal 

We know the location and duration of each disturbance with some certainty. Our search 
algorithm is fed the locations and durations of each disturbance and is then asked to 
remove those disturbances as best as possible. The bulk of our analysis had only one 
disturbance present in the data set. The disturbance function is a quadratic. It easily 
can be shown that the accuracy to which you can fit a quadratic is dependent on the 
number of points available. Therefore the location of the disturbance within our 38-day 
data set should affect our ability to accurately identify the disturbance parameters. 

We examined three disturbance locations in detail. For the rest of this article we 
refer to the three disturbances as gl, g2, and g3. gl has a disturbance start time of 84,000 
seconds (0.97 days) after the beginning of the time series. g2 has a disturbance start time 
of 500,000 seconds (5.8 days) after the beginning of the time series. g3 has a disturbance 
start time of 1,000,000 seconds (11.6 days) after the beginning of the time series. Each 
of these disturbances had the following parameters [p, u, a]i nput = [10~ 17 , 10~ 23 , 10~ 28 ], 
which are referred to as the input disturbance parameters. These values correspond to a 
position change of about 26 nm, a velocity change of about 26 fm/s, and an acceleration 
change of about 0.26 attometer/s 2 of the proof mass. The length of the data interruption 
is ten data points, i.e., 1000 seconds, or about 17 minutes. 

Table E] shows the results of our search algorithm when no signal is present. The 
determination of the position is very good, and the derived velocity and acceleration 
make a smooth residual curve throughout the interruption region. Note that the 
algorithm has found a lower x 2 than obtained by using the input parameters. Figure El 
contains our results for all 35 noise realizations. We can use the standard deviation of the 
derived parameters throughout the 35 noise realizations as an estimate of our error in any 
one derived parameter. The errors are roughly [0.01%, 18%, 0.67%]. Therefore we can 
conclude that our algorithm finds confident values for the position and acceleration, but 
the velocity may not be accurate. Note that in figure El the acceleration is systematically 
higher than the input value, and that the velocity (by correlation is therefore) lower 
than the input value. This slight shift in the acceleration and velocity has allowed our 
algorithm to find a lower \ 2 value. The fact that we have found "incorrect" values 
for the disturbance parameters should not deter us since the uncertainties should not 
depend on the assumed values, and since our main goal is to see if and how this affects 
the science results. 
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Figure 2. Search algorithm results for the three disturbance locations examined, in the 35 realizations of the LISA data stream. The 
dashed lines represent the input parameters. In the case of x 2 , the dashed line is the value of x 2 in the original LISA data stream 
realization. From the collection of 35 values for each parameter we compute the standard deviation. We can take this number as an 
indicator of the error in our derived quantity. This number is printed at the bottom of each plot, for each disturbance location. 
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Figure 3. Search algorithm results for the 35 realizations of the LISA data stream in the presence of a SNR=5 signal only. The dashed 
lines represent the input parameters. In the case of x 2 , the dashed line is the value of x 2 in the original LISA data stream realization. 
The standard deviation of each collection is printed at the bottom of each plot. 
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Table 1. Search results for the three disturbance locations examined with no signal 
present. In all three cases, the search algorithm has minimized \ 2 below that which 
the input parameters give, p* = (p — 10~ 17 ) x 10 20 , u is in multiples of 10 -23 , a is in 
multiples of 10~ 28 . 
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3.3. Signal Identification and Removal 

For our realization of the LISA data stream, an rms strain amplitude of 2 x 10~ 17 
has a signal-to-noise ratio (SNR) of 5 for a 3 (xHz signal. We define the SNR using 
the frequency spectrum of the data set. The SNR is the amplitude of the sine wave 
divided by the root mean squared noise level at the frequency of the sine wave. 
We injected a signal with this amplitude, frequency, and a phase of zero degrees 
into our data stream, so that (A, f, (p) input = (2 x 1CT 17 , 3(xHz, 0°). Without a 
disturbance in the data stream, our search algorithm found these signal parameters 

(A, f, 0) no dis = (2.045 x 10- 17 ,2.996191m.Hz, 4.990°) with Xnodis = 21164. These will 
be the parameters used in our calculation of our supercriterion A 2 . Using the input 
signal parameters we found xLput = 31219. Taking the sum of the squares of the 
difference between the solution and the input signal gives a value of 0.437 x 10 -31 . The 
correlated noise in the data stream "pulls" the solution away from the input signal 
parameters. In effect, our search algorithm has been able to lower x 2 by fitting some of 
the noise. 

Figure El displays the results from our search algorithm for all 35 noise realizations 
in the case of a signal only, with no disturbance. This data shows how well our algorithm 
can find the 3 (xHz signal out of the noisy data. We located the amplitude to within 
1%, the frequency is consistently low by about 2%, and the phase is correspondingly 
high by about 6 degrees. The systematic error in our derived frequency is a result of 
the correlated noise rising as 1/f 3 at lower frequencies. Our algorithm is able to achieve 
a lower x 2 by fitting a slightly lower frequency sine wave. The sine wave frequency is 
anti-correlated with the phase, so the phase is correspondingly high when the frequency 
is low. This systematic error should not deter us in this investigation because we wish 
to find the limitations that disturbances add to our identification of the science signal, 
not necessarily how well we can fit the sinusoid individually. 

4. Effects on Science 

We now examine the full LISA data stream with a monochromatic signal and a 
disturbance. 
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Table 2. Search algorithm results for the three disturbances examined, and 
with a low frequency signal. A is in multiples of 10~ 17 , /* = (/ — 3uHz)/nHz, 
p* = (p — 10~ 17 ) x 10 20 , u is in multiples of 10 -23 , a is in multiples of 10~ 28 , and 
A 2 is in multiples of 10~ 31 . 





a r ct> 


p* u a 


X 2 A 2 


input 
nodis 


2.000 0.000 0.000° 
2.045 -3.809 4.990° 


0.000 1.000 1.000 


31369 0.437 
21164 0.000 


.92 

53 


1.995 -10.754 1.875° 
1.991 -14.507 -0.378° 
1.998 -10.561 -1.322° 


0.190 1.264 0.996 
0.012 1.037 1.009 
-0.170 0.667 1.027 


21449 2.443 
22088 6.196 
22914 4.628 



4-1. Lower Frequency Signal 

A strong low frequency signal can alter the determined disturbance parameters. This 
is why we utilize random restarts of our simplex to prevent becoming caught in local 
minima. In this way we can better determine all six parameters of interest. Table |2] 
contains our search results for the signal and disturbance parameters for the three 
disturbance locations examined. Our search algorithm was again able to minimize y 2 
below that which the input parameters yield. 

Comparing the disturbance values in this table with those in Table ^ we see that 
the derived positions and accelerations are very close, but the velocities are different. 
Comparing the results for all 35 noise realizations we see that in Figure El (which contains 
a disturbance by no signal) there is a systematic low in the velocity, whereas in Figure 0] 
(which contains both a disturbance and a signal) the systematic low has been removed. 
The Fourier transform of a linear slope is a 1/ f power spectrum. A monochromatic sine 
wave will have 1/ / wings. Therefore it is not unreasonable for there to be a correlation 
between the frequency of the sinusoid (which determines the amplitude of the wings) 
and the velocity. This correlation has caused the derived velocities to be different when 
there is a sine wave present. 

Figure 0] contains results for all 35 noise realizations. From these results 
we again compute standard deviations and use these as estimates of the errors 
in each of our derived parameters. The errors in the signal parameters are 
roughly (2.3%, 0.12%, 0.53%). The errors in the disturbance parameters are roughly 
[0.02%, 32%, 1.3%]. With a disturbance present our algorithm has performed worse at 
obtaining the true values of the signal parameters. In addition, the spread in derived 
disturbance parameters has also increased. However, the derived errors from the 35 
realizations imply that the input parameters are just within our solution. This does 
not mean that our algorithm did a good job at locating the signal, rather that the 
disturbance did not significantly push the solution away from the input. 

4-2. Higher Frequency Signal 

Disturbances should have a stronger effect on low frequency sources than high frequency 
sources. At higher frequencies, data disturbances should be less noticeable since there 
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Figure 4. Search algorithm results for the 35 realizations of the LISA data stream with a SNR=5 3 [xRz signal and a disturbance with 
an interruption length of 17 minutes. The dashed lines represent the input parameters. In the case of x 2 , the dashed line is the value of 
X 2 in the original LISA data stream realization. As in the other figures, the standard deviations are printed at the bottom of each plot. 
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are many more cycles present in the data set. Fitting for the disturbance parameters 
with only a high frequency signal present should be much like fitting for the parameters 
with no signal present at all. 

Table 3. Search algorithm results for the three disturbances examined, and with 
a high frequency signal. A is in multiples of 10~ 21 , /* = (/ — 100^Hz)/nHz, 
p* = (p — 10~ 17 ) x 10 20 , u is in multiples of 10~ 23 , a is in multiples of 10~ 28 , and 
A 2 is in multiples of 10~ 39 . 





A f* 4> 


p* u a 


X 2 A 2 


input 
nodis 


1.600 0.000 0.000° 
1.682 -10.050 10.602° 


0.000 1.000 1.000 


31364 1.302 
30692 0.000 


hi 
h2 
h3 


1.687 -9.950 10.413° 
1.689 -10.080 10.778° 
1.684 -10.040 10.586° 


-0.148 0.652 1.011 
-0.047 0.930 1.003 
-0.121 0.750 1.011 


29925 0.366 
30916 0.382 
30185 0.394 



We choose to look at a 100 (J.Hz signal for our higher frequency case. A signal- 
to-noise ratio of 5 implies an rms strain amplitude of 1.6 x 10~ 21 for our realization. 
Our input parameters are now (A, f,(j))- mput = (1.6 x 10~ 21 , 100|xHz, 0°). With 
no disturbance, our search algorithm found the following parameters: (1.682110 x 
10~ 21 , 99.98995n.Hz, 10.60217°) with x 2 = 30692. Table H contains results for this high 
frequency signal in the presence of a disturbance with an interruption of 17 minutes 
in length. Comparing with Table ^ we see that the derived disturbance parameters 
are almost identical. The frequency of the signal has been found to within 10 nHz, or 
0.015%, of the input frequency. The disturbance has had little effect on the ability of 
our algorithm to determine the high frequency signal parameters. 

5. Disturbance Limitations on Operations 

As we have seen in the previous section, isolated disturbances do not significantly affect 
our ability to identify and remove sources from the data stream, given enough data 
before and after the disturbance location. The following question can now be posed: 
How does the frequency of disturbances affect our search algorithm? For this analysis we 
examined only the lower frequency signal case since the identification of higher frequency 
signals is not perturbed much in the presence of a disturbance. 

5.1. Frequency 

Each disturbance contains three parameters to be fitted for. As mentioned earlier, 
the accuracy to which you can determine the parameters depends on the number of 
sample points available. Within a 38-day period, the inclusion of multiple disturbances 
effectively reduces the amount of information available to each disturbance. We can 
simulate the effect of having multiple disturbances by placing one disturbance near the 
end of the data stream. This assumes that all previous disturbances have been fit and 
removed. Table HI shows our results for one disturbance, with L = 1000 seconds, at 
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Table 4. Search results for a low frequency signal with one disturbance. The 
disturbance start time is T before the end of the data stream, and the length L is 
1000 seconds. The last three data sets are g3, g2, and gl from Table [21 A is in 
multiples of 10" 17 , /* = (/- 3|aHz) /nHz, p* = (p - 10~ 17 ) x 10 20 , it is in multiples of 
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4.990° 










21164 


0.000 


1.2 hrs 


9 


uuu 


-0.049 


0.007° 


-0.048 


1.013 


17 


369 


32933 


0.441 


2.4 hrs 


1 


999 


-0.223 


0.029° 


-0.096 


0.849 


19 


351 


30669 


0.454 


6.0 hrs 


1 


998 


-0.666 


0.091° 


-0.241 


0.431 


10 


602 


28215 


0.492 


12 hrs 


1 


996 


-1.237 


0.168° 


-0.311 


0.284 


5 


710 


25756 


0.549 


1.0 day 


1 


989 


-3.552 


0.493° 


-0.145 


0.625 


3 


548 


21438 


0.866 


1.2 days 


1 


991 


-2.790 


0.388° 


-0.404 


0.392 


2 


737 


22334 


0.745 


2.4 days 


1 


985 


-6.858 


1.562° 


0.434 


1.941 


1 


191 


21653 


1.271 


3.5 days 


1 


987 


-4.003 


0.568° 


-0.371 


0.194 


1 


457 


24205 


0.942 


5.0 days 


1 


988 


-5.911 


2.318° 


-0.131 


0.612 


1 


194 


21772 


0.747 


7.0 days 


2 


001 


-3.662 


1.557° 


-0.323 


0.166 


1 


171 


24628 


0.519 


9.5 days 


1 


986 


-6.369 


2.492° 


-0.271 


0.547 


1 


086 


21947 


0.796 


14 days 


1 


929 


-8.449 


1.817° 


0.461 


1.699 





963 


25784 


1.905 


19 days 


2 


021 


-6.491 


2.646° 


0.302 


1.535 





976 


22295 


0.695 


23 days 


1 


978 


-8.641 


3.191° 


0.173 


1.395 





987 


22349 


1.132 


26 days 


1 


998 - 


-10.561 


-1.322° 


-0.170 


0.667 


1 


027 


22914 


4.628 


32 days 


1 


991 - 


-14.507 


-0.378° 


0.012 


1.037 


1 


009 


22088 


6.196 


37 days 


1 


995 - 


-10.754 


1.875° 


0.190 


1.264 





996 


21449 


2.443 



various distances from the end of the data stream. With a small number of points 
after the data interruption it is hard to determine the acceleration parameter. This 
implies that the assumption that we removed all other disturbances exactly is wrong — 
we probably would not have been able to determine the acceleration well. Therefore we 
must fit for all of the disturbances at once since the acceleration parameter will strongly 
affect the data at large times. 

As a preliminary investigation we examined several cases with multiple 
disturbances. We simulated data streams with the maximal number of disturbances 
possible all equally spaced by T, and of length L = 1000 seconds. Our results of this 
analysis are in Table EJ Note that our search algorithm was not designed to handle more 
than 20 disturbances, therefore instead of having 31 equally spaced disturbances in the 
T = 1.2 days data set and 37 disturbances in the T = 1 case, we only have 20. 

It appears that our algorithm has trouble with disturbance frequencies of more 
than one every 3.5 days. The algorithm not only has trouble fitting for the frequency 
and phase, but the velocities and accelerations for most of the the disturbances are 
much worse — this is evident by the large values of \ 2 ■ Some refinement of the search 
algorithm for multiple gaps may help to fit the velocities and accelerations better for a 
high frequency of disturbances. 
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Table 5. Search algorithm results in the presence of multiple disturbances. # is the 
number of disturbances present in the data stream, separated in time by T, and all 
of length L — 1000 seconds. Notice the degradation in the derived signal parameters 
and the value of x 2 f° r disturbance frequencies more than one every 3.5 days. A is in 
multiples of 10~ 17 , /* = (/ - 3uHz)/nHz, and A 2 is in multiples of 10~ 31 . 



T 


# 


A 


/* 


<f> 


x 2 


A 2 


input 





2.000 


0.000 


0.000° 


31218 


0.437 


nodis 





2.045 


-3.809 


4.990° 


21164 


0.000 


1.0 day 


20 


2.447 


-403.809 


12.292° 


4991500 


287.479 


1.2 days 


20 


2.093 


-192.644 


12.261° 


3853900 


311.622 


2.4 days 


15 


1.991 


-2.163 


-1.474° 


819324 


1.329 


3.5 days 


10 


2.046 


-3.735 


5.180° 


29346 


0.002 


5.0 days 


7 


1.993 


-6.640 


3.392° 


24459 


0.561 


7.0 days 


5 


2.044 


-5.313 


4.941° 


23877 


0.048 


9.5 days 


3 


2.076 


3.460 


3.027° 


23818 


0.277 


14 days 


2 


2.018 


-3.579 


0.781° 


21993 


0.708 


19 days 


1 


2.108 


-1.311 


2.381° 


21393 


0.222 



6. Conclusions 

We have presented an algorithm which identifies monochromatic signals in the presence 
of data disturbances in a LISA-like data stream. The results here are encouraging on 
the ability of removing the disturbances without affecting the science results. However, 
we found complications with data sets containing disturbances which are more frequent 
than about one every 3.5 days. There are a number of simplifications in our data 
stream, for instance the use of monochromatic signals and instantaneous drop out and 
recovery on the edges of interruptions. As mentioned in ^3.21 longer data sets will allow 
all parameters to be fit more accurately. Therefore the results from our short 38 day 
long data set analysis should be taken as a guide for longer data sets. 
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